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Dielectric Correction to the Chiral Magnetic Effect 

Kenji Fukushima and Marco Ruggicri 
Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

We derive an electric current density j em in the presence of a magnetic field B and a chiral 
chemical potential /15. We show that j has not only the anomaly-induced term oc [i^B (i.e. 
Chiral Magnetic Effect) but also a non-anomalous correction which comes from interaction effects 
and expressed in terms of the susceptibility. We find the correction characteristically dependent on 
the number of quark flavors. The numerically estimated correction turns out to be a minor effect 
on heavy-ion collisions but can be tested by the lattice QCD simulation. 



PACS numbers: 12.38.Aw,12.38.Mh 
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The Chiral Magnetic Effect (CME) 0, Q is a general 
mechanism to induce an electric current j em along the 
direction of an external magnetic field B for systems that 
have a non- vanishing chirality charge N$ = Nr — Nr,. 
Using the chiral chemical potential /is to express N$ ^ 
in the grand canonical ensemble one can find the electric 
current density given simply as [H 



e 2 /i 5 -B 
2tt 2 



(1) 



which originates from the quantum anomaly and thus is 
an exact relation insensitive to any infrared scales such 
as the particle mass, temperature, etc. This type of 
anomaly relation is quite generic and also discussed in 
various contexts Particularly, in prior to the CME, 

a dual situation had attracted attention; it is the axial 
current j 5 that is induced by B if the quark chemical 
potential \i q is non-zero, which is parallel to each other 
under the replacement /Ug -H- \i q and j -f-> j 5 . The re- 
sulting axial current may be realistic in the cores of com- 
pact stellar objects [H or in some condensed matter sys- 
tems @ . Generally speaking, the structure of Eq. (JXJ) can 
also arise from low-energy effective descriptions. That is, 
from symmetry reason, the anomalous current is to be ex- 
pressed as jg m = —{d y (f)F tlL ' where <j> is a pseudo-scalar 
field 0] • Under the magnetic field B l = F° l the cur- 
rent is thus proportional to B and 8q4> that translates 
into ^5 of Eq. (JTJ). This gives a physical interpretation 
of /Ltg as a time derivative of the pseudo-scalar conden- 
sate. Equivalently one can also say that /X5 appears as a 
result of a time derivative on the strong 0-angle parame- 
ter Q which leads to pseudo-scalar meson enhancement 
through mixture of the a and 770 condensates 0, Q . 

It is tremendously important to verify Eq. ([T]) in or- 
der to quantify the CME in such a way that theoretical 
predictions can be compared to experimental measure- 
ments [f]| . At the same time the estimate of back- 
ground (non-topological) effects 12 -15| is indispensable 
for experimental confirmation of the CME. The purpose 
of this Letter is to point out that the induced current ([T]) 
itself must receive a non-anomalous correction from in- 
teraction effects even at the level of the mean-field quasi- 
particle description. The notable feature of the correc- 



tion is that the coefficient in Eq. (p}, which seems to be 
protected by anomaly, is modified by a factor. Existence 
of such correction is a novel insight in theory. 

In this Letter we do not directly consider the interac- 
tion mediated by gauge bosons but instead make use of 
an effective form of the interaction in terms of fermionic 
degrees of freedom. It should be a reasonable approxi- 
mation to utilize the current-current interaction, 



(2) 



as long as the typical interaction energy is lower than the 
gauge boson mass, which is the case for Fermi's effective 
theory of the weak interactions. In QCD the one-gluon 
exchange can result in an interaction form of the color 
current in which the color generators are inserted. It 
is then possible to extract a specific form ^ out from 
the Fierz transformation. Accordingly Gy should be of 
order ~ g /Mi where g is the gauge coupling constant 
and M g is the gauge boson mass. [Gluons should be 
massive non-perturbatively though their mass is zero at 
the QCD Lagrangian level.] For the moment we treat 
Gy as just a parameter and will plug a concrete value in 
later discussions. 

The vector interaction ([2]) is special in QCD; it is in- 
variant under chiral rotations and thus its presence is 
naturally anticipated. Introducing a mean-field j M = 
(•07^) the interaction © is decomposed into £y — > 
—Gvj z2 + '2Gvj z ipJ 3 ''P, then, apart from quantum gauge 
fluctuations, we have the kinetic term; £km = ^{il^^ 1 — 
7\/ + /i57°7 5 )'0, where M is a mass which may be dynam- 
ically generated, but for the present purpose the micro- 
scopic origin of M is irrelevant. In writing the above 
we have introduced an effective (classical) gauge field in 
£>M = - iA* given by 



A x = A v = , A z 



-2G v j> 



(3) 



for which the coupling constant is chosen as unity. [This 
is a convenient choice for later generalization to QCD 
with quark flavors having different electric charges.] The 
grand potential (divided by the volume V) can be read 
from the zero-point oscillation energy in addition to the 
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mean-field condensation energy, i.e. 

n/v = Gyf 2 -&Y, a ^f T 1 ^) , (4) 



s . k 



where s is the spin and k refers to the Landau level. The 
spin degeneracy factor is taken care of a s k defined as 



J s,+1 



a s k 



for 
for 
for 



k = 0, 
k = 0, 



eB > 
eB < 



(5) 



Now the quasi-particle dispersion relations are derived 
from the eigenvalues of the Dirac operator which are Q ; 



uj 2 s = M 2 + [|p|+sgn(p 2 )s / u 5 ] 2 , where \p\ 2 = (p z +A z ) 2 + 
2\eB\k with k being a non-negative integer to label the 
Landau level. 

The stationary condition for Q with respect to j, i.e. 
dfl/dj = 0, yields a self-consistent condition to deter- 
mine the current density, 







d(Q/V) 
OA 



(6) 



We note that the derivative with respect to A acts only 
on the latter term of (HI that explicitly depends on A. 

The expression ^ looks like a standard one utilized 
in Ref. [2| but the essential difference is that the deriva- 
tive should be evaluated at non-zero A which is specified 
by (0 and thus j should be solved in a self-consistent 
manner. 

By taking the derivative explicitly on Eq. (j4} only the 
surface terms contribute to the current. Here we assume 
some regularization scheme and introduce a ultraviolet 
momentum scale in such a way that A 2 +2\eB\k = A 2 . 
Then the surface terms are picked up as 



d(Sl/V) 



dA z 



\eB\ 
An 2 



fe=l s 

with the energies at p z = ±Afc, namely, 



(7) 



9± = VM 2 + (|±Ao + ^|±sgn(eS)/i 5 ) , 

«A„±/± sgn(eB) M5 , (8) 

for the Landau zero-mode k = and the second line is an 
approximation valid for sufficiently large A. In the same 
way we have 



f s ± = ^M 2 + (v/(±A fc +^) 2 +2| eJ B|fc + sus) 

« A ± ^A z 
A ' 



(9) 



for the Landau non-zero modes k > 0. From our defini- 
tion of Afe and A it is obvious that fc max = [A 2 /(2|ei?|)J . 

It should be noted that the first contribution in ([7]) 
involving g± leads to the electric current (fTJ) in the limit 



of vanishing A. The latter term involving f s ± is simply 
zero if A is absent. 

When A is sufficiently larger than \eB\, the sum over 
k can be well approximated as J^fcO^-fc/A-) « A 2 /(3|e£?|). 
Using this we reach, 
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|efl| 
2tt 2 



sgn(e5)/i 5 - 2Gv[ 1 



2A 2 
3|eS| 
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(10) 



From the above one might be able to solve j, though 
the divergent term still remains. The expression is not 
physically meaningful yet as it is. One needs to formulate 
the renormalization procedure that we address in what 
follows below. 

Here it should be mentioned that our result (TTU)) is 
quite analogous to what is discussed in Ref. 0] . The for- 
mulation may look different since no current-current in- 
teraction was introduced in Ref. Q but familiar Nambu- 
Jona-Lasinio (NJL) type (scalar and pseudo-scalar) four- 
fcrmion interactions were considered with the Dyson- 
Schwingcr equation which goes beyond the present mean- 
field analysis. Nevertheless our treatment suffices to 
grasp the essential point of Ref. 0. That is, the in- 
teraction term like Eq. © would be always induced by 
mixing interactions from scalar and pseudo-scalar chan- 
nels, which might be the case in the Dyson-Schwinger 
calculation. 

For the purpose of renormalization let us consider the 
following susceptibility defined by 



C 



d 2 {Sl/V) 



dA z 



\eB\ 
2^" 



1 



2A 2 
3|efl| 



(11) 



which is still divergent. Because this susceptibility can be 
interpreted as the gauge boson (screening) mass, it must 
become vanishing in the vacuum (i.e. for B = 0) so that 
the gauge invariance is maintained. This imposes the fol- 
lowing requirement; C — > for B — > 0, which sets a nat- 
ural condition for minimal subtraction. In the definition 
of C in Eq. (fTTj) the second term is divergent as A 2 /(37r 2 ), 
which must be subtracted through renormalization. Con- 
sequently the renormalized susceptibility should be given 
by the finite first term; 



Cr 



\eB\ 
2^" 



(12) 



This result is consistent with the conclusion in Ref. [13| 
in which an independent derivation of the susceptibility 
(|12[) is given by means of the linear response formula with 
the current generation due to anomaly [l6| . Since the 
derivation in Ref. [L3| is free from ultraviolet divergence, 
this confirms the validity of Eq. (TT2"|) . 

It is straightforward to solve (ITU1) with respect to j z 
with C replaced by Cr. Wc finally find, 



1 



1 + 2G v Cr ' 2ir 2 



(13) 
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from which, apart from flavor complication in QCD 
which will be discussed later, the electric current density 
is simply j cm = ej. If we interpret the coefficient in front 
of B as the chiral magnetic conductivity [l7[, Eq. (fTB")) 
means that the conductivity should be corrected by the 
in-medium (i.e. magnetic field induced) (generalized) di- 
electric constant k [18[. We here define K by 



7.00 



«-icm(Gv =0) 



(14) 



We note that k deduced from Eq. (fl"3"j) is a standard for- 
mula representing the screening effects by vacuum polar- 
ization [18|. 

As wc pointed out, the current-current term like 
Eq. ([2]) should be generally present as a result of non- 
perturbative interactions, and thus Eq. ([TJ is no longer 
the exact answer in the fully interacting case. In this 
sense we would claim that the holographic calculations 
as in Ref. [lj| may miss some back-reaction because the 
correction factor ([T3]) has not been found there. This dif- 
ference might explain subtleties on the holographic chiral 
magnetic current which are still under dispute 111 20j |. 

Let us generalize our result (fl~3j) to the QCD prob- 



lem. Then, the complication comes from flavor degrees of 
freedom. In what follows we will address the one-flavor, 
two-flavor, and three-flavor cases in order. 

One-flavor case: For the case with only one flavor the 
modification is only the color factor N c = 3 and also that 
e may be replaced by the quark electric charge q. Then, 
the dielectric constant is almost trivially, 



K= 1 



— G v \qB\ 

7T Z 



(15) 



which is unity at Gy = and goes to zero as Gy — > oo. 

Two-flavor case: The situation is drastically changed 
if there arc multiple flavors with different electric charges. 
Now we consider a system having u quarks with the elec- 
tric charge q u = (2/3)e and d quarks with qd = (— l/3)e. 
The susceptibility becomes flavor dependent as 



C 3 



Nc\q f B\ 
2ir 2 



(16) 



with which the current is expressed in the same way as 
Eq. (|T3J) as follows; 



1 



Ms qfB 

1 + 2G V J24 " 2 ^ 



(17) 



'R f 
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Up to this point the generalization is just straightfor- 
ward. In the presence of non-degenerate flavors the elec- 
tric current is no longer proportional to j, but one has to 
compute j cm as the flavor summation of qj times current 
contribution from each flavor sector. That is, the electric 
current density should be read from 



= N, 



E 



2tt 2 



2G v jY,QfC 



f 

R ■ 



(18) 




FIG. 1. The dielectric correction coefficient as a function of 
GylcBI in the two-flavor case. 



After some calculations wc find, 



K = 1 



12 

5^ 



Gy\eB\ l + —G v \eB\ 



(19) 



This result (|T9j) behaves different qualitatively from 
Eq. (fTB"]) . One can easily see that k asymptotically 
approaches a finite number 4/5 = 0.8 in the limit of 
Gy — > oo and k never goes to zero. To show this clearly, 
we make a plot for the above function of GV|e-B| in 

Fig. [TJ It is obvious from the figure that k slowly decays 
to the asymptotic value 0.8 and thus the dielectric cor- 
rection is only a minor effect in contrast to the one-flavor 
situation. 

Three-flavor case: It is interesting to think of the 
ideal case with three flavors. Because the chiral magnetic 
current has the anomaly origin and is independent of the 
quark masses, the three-flavor case might be realistic if 
/X5 is large enough. Let us imagine that the system has s 
quarks with q s = (— l/3)e in addition to u and d quarks. 
Then the situation is totally changed again. In view of 
Eq. p7|) . j is proportional to q/ which becomes van- 
ishing for the three-flavor case; q u +qd + qs = 0. Because 
j is zero, there is no correction appearing at all, and thus 
identically K = 1. This result is intuitively understand- 
able if quark masses are degenerate; the system is then 
automatically electric-charge neutral. Hence, there is no 
coupling between the baryon current j and the electric 
current j em and thus no back-reaction from their entan- 
glement. It is, however, non-trivial that the conclusion 
of no correction for the three-flavor case holds regardless 
of whether quark masses are degenerate or not. 

Now let us plug concrete numbers in our final expres- 
sion to see how large/small the correction is in specific ex- 
amples. We first need to determine a value of Gy. From 
the discussions below ([2]) we can postulate Gy ~ g 2 /Mg 
and let us choose g = 2 (a s ~ 0.3) and M g = 0.8 GeV 
(that is roughly a half of the glueball mass) here. Then we 
have a rough estimate as Gy ~ 6.3 GeV~ 2 . We can make 
it sure that this is a reasonable estimate from the empir- 
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ically adopted Gy in the NJL model; Gy = 0.2 ~ 0.5Gs 
where Gs is the four-fermion coupling in the scalar and 
pseudo-scalar channel and fixed as Gs — 9.2 GcV~ 2 to 
reproduce the pion mass and decay constant (2lj . This 
is not far from our estimate Gy ~ 6.3 GeV -2 . 

In the heavy-ion collision it would be convenient to 
express the magnetic field strength \eB\ in the unit of 
the pion mass squared m 2 instead of gauss. From the 
UrQMD simulation |ei?| is evaluated as a few times 



for the RHIC energies [22j. If we use the value \eB\ = m 2 , 
then, we find Gy\eB\ = 0.11. This is a small number and 
the dielectric constant stays close to unity for any case of 
flavor number. Therefore, fortunately, we can conclude 
that the dielectric correction from back-reaction is only 
minor and practically negligible for phenomenology. 

Finally let us make a comment on a possible appli- 
cation of our result to the lattice-QCD simulation. The 
chiral magnetic effect has been investigated in the lattice- 
QCD simulation 23, 24] with extremely strong magnetic 
fields. For example, in Ref. |23j], the applied magnetic 
field can be as strong as \eB\ ~ GeV 2 , which is of order 
hundred in the unit of m 2 . If we use Gy ~ 6.3 GeV -2 
then Gy|eB| ~ 6.3 for \eB\ ~ 1 GeV 2 . According to our 
expressions the one-flavor system would lead to a sub- 
stantial suppression factor k ~ 0.34. Even in the two- 
flavor case the suppression is a sizable effect; k ~ 0.87. 
Therefore, it should be possible to confirm the existence 
of such dielectric corrections as discussed here using the 
lattice-QCD simulation. 

The lattice-QCD simulation opens an intriguing pos- 
sibility that Gy may be determined from k. In fact the 
determination of Gy provides us with very useful infor- 
mation on the QCD phase diagram. Especially it cru- 
cially depends on Gy whether the chiral phase transition 
can become of first order at finite density and whether 
the QCD critical point can exist on the phase diagram. 
Once a finite baryon chemical potential is turned on, of 
course, the sign problem hinders the simulation. Never- 
theless the two-color two-flavor simulation is still feasible 
even at finite density, which may give Gy as a function 
of density, if the precise determination of chiral magnetic 
current is possible from the lattice data. 

In summary we computed the back-reaction coming 
from the vector interaction, which should result in a di- 
electric correction on the chiral magnetic current even at 
the mean-field level. Our final expressions show that the 
qualitative behavior of the correction strongly depends 
on the relevant number of flavors in the system. The 
one-flavor case has a substantial suppression on the chi- 
ral magnetic current due to screening effects, while the 
two-flavor case has only a minor modification however 
strong the vector interaction is. There is no correction at 
all for the three-flavor case. It should be possible to quan- 
tify the correction in the lattice-QCD simulation, which 
in turn would give useful information on the strength of 
the effective vector interaction. 
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